Nuclear multifragmentation, percolation and the Fisher Droplet Model: common 

features of reducibility and thermal scaling 
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It is shown that the Fisher Droplet Model (FDM), per- 
colation and nuclear multifragmentation share the common 
features of reducibility (stochasticity in multiplicity distribu- 
tions) and thermal scaling (one-fragment production proba- 
bilities are Boltzmann factors). Barriers obtained, for cluster 
production on percolation lattices, from the Boltzmann fac- 
tors show a power-law dependence on cluster size with an 
exponent of 0.42 ± 0.02. The EOS An multifragmentation 
data yield barriers with a power-law exponent of 0.68 ± 0.03. 
Values of the surface energy coefficient of a low density nu- 
clear system are also extracted. 



Since the earliest observations of nuclear multifrag- 
mentation (the break up of excited nuclei), the Fisher 
Droplet Model (FDM) @] and percolation models ||] 
have been employed in attempts to understand this phe- 
nomenon. The FDM enjoyed early success in predicting 
power-law distributions in fragment masses at the critical 
point in a liquid- vapor diagram jsj. Percolation models 
also predicted a power-law distribution in fragment sizes 
near the critical point Q, Both models still enjoy 
great popularity and have been employed in the analy- 
sis of Au multifragmentation data obtained by the EOS 
Collaboration §,§,§, 1], 0. 

Other analyses of multifragmentation data have shown 
two empirical properties of the fragment multiplicities 
which have been named reducibility and thermal scaling 
0, H, [Hi, Reducibility refers to the observation 
that for each energy bin, E, the fragment multiplicities, 
N, are distributed according to a binomial or Poissonian 
law. As such, their multiplicity distributions, P^, can 
be reduced to a one-fragment production probability p, 



according to the binomial or Poissonian law: 
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where M is the total number of trials in the binomial dis- 
tribution. The experimental observation that Pn could 
be constructed in terms of p was considered evidence for 
stochastic fragment production, i.e. fragments are pro- 
duced independently of each other. Experimental frag- 
ment multiplicity distributions were observed to change 
from binomial to Poissonian under a redefinition of frag- 
ment from 3 < Z < 20 to individual charges, Z [ p^ . 

Thermal scaling refers to the feature that p be- 
haves with temperature T as a Boltzmann factor: p cx 
exp{—B/T). Thus a plot of Inp vs. 1/T, an Arrhcnius 
plot, should be linear if p is a Boltzmann factor. The 
slope B is the one-fragment production barrier. Analy- 
ses of multifragmentation distributions along these lines 
have demonstrated the presence of these features 0] and 
have led to the extraction of barriers [0 . Controversy 
has surrounded this type of analysis regarding both the 
physical existence of these features and their significance, 
mostly within the framework of dynamical vs. statistical 
origins of multifragmentation |13|. 

In this work several important points will be made: 
The FDM inherently contains reducibility and thermal 
scaling. Since percolation reduces to the FDM, it exhibits 
reducibility and thermal scaling. Thus percolation pro- 
vides a simple mathematical model that fully manifests 
these two features. Arrhcnius plots for percolation can 
be used to extract barriers. The barriers have a power- 
law dependence on cluster size. Analysis of the EOS Au 
multifragmentation data verifies reducibility and thermal 
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scaling. The extracted barriers also obey a power-law de- 
pendence on fragment mass. 

The FDM and its forerunners jl^] , are based on the 
equilibrium description of physical clusters or droplets. 
The mean number of droplets of size A was written as: 
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where A/i — fJ- — fJ.i and /i and fii are the actual and liquid 
chemical potentials respectively. For fJ. < m (gas), {Na) 
falls to zero with increasing A. For /i > /i/ (liquid), (Na) 
increases with A. To better describe the distribution for 
intermediate values of A, Eq. (||) was modified to include 
the surface of the droplets: 



(Na) c>c exp 



AAfi c{T)A 



2/3 



T 



T 



(3) 



where c(T) is the surface free-energy density. For ^ < jii, 
{Na) falls to zero with increasing A. For n > ni, the 
terms in the exponential compete, leading to an early 
decrease in {Na) with A, followed by an increase. 

To account for the properties near criticality. Fisher 
introduced an explicit expression for the surface free en- 
ergy and a topological factor resulting in an expression 
for the normalized droplet distribution: 
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where: Aq is the size of the system; is a normalization 
constant depending only on the value of r ]l7[ |; r, the 
topological critical exponent, depends on the dimension- 
ality of the system with origins that lie in considerations 
of a three dimensional random walk of a surface closing 
on itself, for three dimensions 2 < r < 3; cqeA'^ is the 
surface free energy of a droplet of size A; cq is the surface 
energy coefficient; a is the critical exponent related to the 
ratio of the dimensionality of the surface to that of the 
volume; and e = (Tc — T)/Tc is the control parameter, a 
measure of the distance from the critical point, Tc- 

From this outline it is apparent that the FDM exhibits 
the features of reducibility and thermal scaling. The dis- 
tribution in droplet size is Poissonian by construction: in 
the FDM each component of droplet size A is an ideal gas 
without the canonical constraint of overall constituent 
number conservation. The resulting grand canonical dis- 
tribution is Poissonian. Thus, cr^ = {Na), *-e- Poisso- 
nian reducibility. 

Thermal scaling is obvious in the FDM when Eq. (^) 
is written as follows: 
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It is clear that linearity with l/T (thermal scaling in 
an Arrhenius plot) extends to and beyond the critical 



point, and the slope of the Arrhenius plot gives the T = 
surface energy coefficient of the droplet. 

Percolation models are characterized by a constant en- 
ergy per bond. The bond- breaking probability, pbreak, 
is amenable to a straightforward statistical mechanics 
treatment. Such a treatment reveals that in the limit 
of T — > oo, Pbreak — > 1/2, indicating that the range of 
Pbreak covered by Eq. (|^) is half the usual range discussed 
in percolation theory; < pbreak < 1- Such a literal 
thermodynamical treatment therefore excludes the criti- 
cal point, pc, of many types of percolation systems from 
thermodynamic consideration. 

However, percolation phenomena, with a geometrical 
phase transition, share with thermal critical phenomena 
the important features of scaling, universality and renor- 
malization group as well as other deep connections [ p^ . 
For example, the scaling behavior observed in percola- 
tion clusters can be described by the FDM when pbreak 
replaces T in Eq. (||) and the control parameter becomes 

e = {pc - Pbreak) /Pc ^■ 
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FIG. 1. The scaled cluster (fragment) distribution plotted 
as a function of the scaled control parameter for clusters (frag- 
ments) of size (mass) A. The solid line shows a fit to the FDM 
for (a) percolation ((b) Au multifragmentation). 

To demonstrate the scaling of percolation clusters a 



plot is made of the scaled cluster distribution, n 



scaled 



/qqA as a function of the scaled control parain- 

^scaled 



eter, e^'='"'=° = A-'e/pbreak- See the left panel of Fig. ||. 
Data over a wide range in A and e are seen to collapse. 
Cluster distributions used in this analysis were generated 
on a simple cubic lattice of side six . 

Fitting 71^^°'^'' as a function of ^ '^'^ for e > and 
leaving cq and exp [AA^/pbreak] as free parameters in 
Eq. (|) gives cq = 2.34 ± 0.03 and exp [AA^i/pbreak] = 
0.95 ± 0.01, i.e. the bulk factor is unity. At the crit- 
ical point, e = 0, the collapsed distribution takes the 
value of one indicating that the cluster distribution fol- 
lows a power law. Away from the critical point, the clus- 
ter distribution predominantly follows the surface term 
in Eq. (^). The FDM does not describe the behavior of 
clusters for e <C 0. Other forms for the EDM's surface 
factor have been suggested to describe cluster behavior 
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on both sides of the critical point |^ , . 

Also shown in Fig. l| is a plot of the EOS Au multifrag- 
mentation d ata fl, [j7|, ||], [|lO|. Here the substitution 
of \feF = ^ E* /Aq for T has been made resulting in a 
control parameter of e = (\/eJ — which for a 

degenerate Fermi gas reduces to (Tc — T)/Tc- The exci- 
tation energy normalized to the mass of the fragmenting 
remnant, e* in MeV/nucleon, excludes collective effects 
1^. The location of the critical point, e*, and values of 
the critical exponents, a and t, were determined previ- 
ously p, §, H, [§, 0. Fitting nY^""^ as a function 
of e^'^'^ for e > and leaving cq and exp [j4A/j,/\/e^] 
as free parameters in Eq. (||) gives cq = 6.4 ± 0.6 MeV 
(via E* = aT^ with a = Aq/U) and exp [AA^J,/^/e*] ^ 
0.8 ±0.1, i.e. the bulk term is consistent with A/i w 0. 
The surface energy coefficient cq is of a somewhat dif- 
ferent nature than the semiempirical mass formula pa- 
rameter {as ^ 17 MeV for T = Q, p = po) or estimates 
for low density nuclear systems (os ~ 6 MeV for T ~ 3 
MeV, p ~ po/3) [|o). The coefficient cq is temperature 
independent; the temperature dependence is given as cge. 

Since the FDM has been shown to contain the features 
of reducibility and thermal scaling and since the scaling 
inherent in the FDM describes percolation, it should be 
possible to observe reducibility and thermal scaling in 
percolation cluster distributions. 
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FIG. 2. Ratio of the variance to the mean number of clus- 
ters (fragments) of size (mass) A as a function oi pbreak (e*) 
for percolation (a) (Au multifragmentation (b)). Error bars 
are statistical. Vertical dashed lines show the location of the 
critical point. 

To address the question of reducibility in percolation, 
cluster multiplicity distributions for bins in pbreak are 
considered. The ratio of the variance to the mean, 
a\/ (Na), of the multiplicity distribution for each cluster 
of size A is an indicator of the nature of the distribution. 
In Fig. ^, such a ratio is shown as a function of ptreak- 
The observed ratio is near one (Poissonian limit) over 
the range oipbreak- Within experimental errors, similar 
behavior is observed for the Au multifragmentation data. 

Examples of multiplicity distributions with Poissonian 
curves calculated from percolation (iV^) are shown in the 



left panel of Fig. g. Poissonian distributions reproduce 
percolation cluster distributions over two or three orders 
of magnitude for all A values; Poissonian reducibility is 
present in percolation. 

Fig. H also shows the multiplicity distributions for 
Au multifragmentation compared to the calculated Pois- 
sonian curves. The agreement between the measured 
and computed distributions confirms the presence of re- 
ducibility in the Au multifragmentation data. 
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FIG. 3. Multiplicity distributions, Pjv, for clusters (frag- 
ments) of size (mass) A as a function of for bins in 
Pbreak (e*) for pcrcolation (Au multifragmentation) on the 
left (right). The lines are Poissonian distributions calculated 
from the measured (Na). 

To verify thermal scaling in percolation, the average 
yield of clusters of size A and its dependence on Pbreak are 
considered. The presence of thermal scaling should mani- 
fest itself through a Boltzmann factor. Again, the substi- 
tution Pbreak for T IS made in accordance with standard 
percolation theory and In (n^) is plotted as a function 
of 1/ Pbreak (Arrheuius plot). See Fig. ^. In most cases 
the Arrhenius plots for individual clusters of size A are 
linear over two orders of magnitude. Thus, thermal scal- 
ing is verified for percolation. The observations that re- 
ducibility and thermal scaling are already present in such 
a simple model suggest that they are deeply rooted fun- 
damental features of multifragmentation processes rather 
than being epiphenomena of complex systems. 

The Boltzmann factor indicates that the slope of an 
Arrhenius plot represents the barrier B associated with 
the production of a cluster. For an interpretation of B in 
percolation, the Boltzmann factor is equated with Eq. (^) 
yielding a power law relating B to the size of a cluster: 
B = cqA"' when A/i « 0. Fitting the extracted barriers 
B (slopes in Fig. H) as a function of A gives an exponent 
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FIG. 4. Normalized average cluster (fragment) multiplicity 
as a function of 1/pbreak (l/\/e^) for clusters (fragments) of 
different size (mass) , A, for percolation (a) (An multifragmen- 
tation (b)). Solid lines show Arrhenius fits. Vertical dashed 
lines show the location of the critical point. 

equal to 0.42 ±0.02 in agreement with the accepted value 
of (7 = 0.45 for 3D percolation. See Fig. ||. The constant 
of proportionahty of the power law gives another measure 
of the surface energy coefficient cq = 2.42 ± 0.03. 
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FIG. 5. The power law relationship between the Arrhenius 
barrier, B, and the cluster size (fragment mass) A for perco- 
lation (a) (multifragmentation (b)). The solid line is a fit to 
the points with A > 1. 

When the Arrhenius analysis is performed on the Au 
multifragmentation data, the results are quaUtatively 
similar, but quantitatively distinct. Arrhenius fits of 
In(n^), now plotted against are linear over an 

order of magnitude or more. See Fig. m The bar riers 
extracted here can be converted into units of MeV via 
the energy-temperature relation for a degenerate Fermi 
gas. Following the same analysis as for percolation, a 
fit was made of B vs. A (see Fig. ||) which yielded 
a — 0.68 ± 0.03, in agreement with previously deter- 
mined EOS Au multifragmentation values [Q, |jl^, ||2^] 
and Co = 6.8 ± 0.5 MeV. 



In summary, the above effort illustrates: 

• the presence of reducibility and thermal scaling in 
the FDM, percolation and the EOS Au multifrag- 
mentation data (the latter two shown empirically) ; 

• the relationship between the FDM and percolation 
via a scaling analysis that also yields an estimate 
of the surface energy coefficient; 

• that the barriers obtained from percolation follow a 
power-law dependence on cluster size with an expo- 
nent that agrees with the accepted 3D percolation 
value and gives another (consistent) estimate of the 
surface energy coefficient; 

• the collapse of the EOS Au fragment distributions 
in accordance with the FDM yielding an estimate 
of the surface energy coefficient of a low density 
nuclear system; 

• that the barriers obtained from EOS Au multifrag- 
mentation data follow a power-law dependence on 
fragment mass with an exponent near the expected 
value (2/3) and close to the 3D Ising universality 
class value (0.64) and gives another (consistent) es- 
timate of the surface energy coefficient. 
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